X 



PROCEEDINGS 

°F SCIENCE 



o 
o 

(N 

> 

o 
O 
in 



> 
oo 
(N 
O 
O 

in 
o 

> 

Oh: 



Towards an algebraic approach to the discretisation 
of fermions. 



Vivien de Beauce*^ 

Department of Physics, Hokkaido University, Sapporo, 060-0810, Japan 
E-mail: |debeauce@particle . sci . hokudai .ac.jp 



A discretisation scheme for differential geometry is applied to the problem of constructing lattice 
actions, the naive and staggered action are thus derived. It is found that after specifying an ansatz 
for the space of fields, the corresponding lattice action is obtained. The gauging procedure, and 
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on the QED axial anomaly are made, for the theory in which the lattice projection operator is not 
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Recently, renewed interest has been directed to model building of lattice theories in relation 
to chiral symmetry The Ginsparg-Wilson relations have provided an optimal way to break 
chirality (P. Hasenfratz), even preserving a remnant of the symmetry (Liischer). The axial anomaly 
in QED is reproduced coiTcctly. Also, the Poincare lemma on the lattice may be written, thus, 
encoding algebraic topology in the formulation. Such considerations suggest the need for a more 
algebraic formulation of lattice (although not using the non-commutative geometry, see the discus- 
sion in [Q]). We argue that differential geometry should also be captured. 

Our starting point is a 2D squai^e complex K (the treatment in 4D is almost identical), together 
with the chain space, i.e linear combinations of vertices, edges, squares (the setup of homology). 
We start with the Dirac-Kahler (DK) equation which is differential geometric. The task is then to 
discretize the various operations required, which are d, the exterior derivative, A the wedge product, 
-k the Hodge star operator, and the interior product All these operations play an important part 
in the model. The first lattice transcription of the DK equations, based on chains was done by 
Becher and Joos A central feature of that method, is that although the fermion field is a chain, 
the coefficients are placed at vertices only (i.e on the lattice), requiring the introduction of lattice 
translation operators in order to satisfy some of the defining algebraic relations of the operators 
(a, d, *). The Dirac equations arising from the DK equation could not be decoupled, and so one 
was left with the lattice DK fermion with its multiplicity of "flavors". Recently it was recognized 
that the Whitney map W of algebraic topology makes for a neater discretization in which there is 
no need for displacement operators [Q], but the underlying lattice plays no special role, so a direct 
comparison with lattice actions was not obvious. In the meantime, a solution to the problem of 
defining a satisfactory discrete interior product was given [||], the key step is the introduction of a 
product X on the space of chains, so that topology is captured in the space of chains as usual, while 
an element of continuous geometry is captured in the product space, in a consistent way. 

In this talk' we revisit the problem of fermions with this new technique in mind. A new discrete 
Hodge star will also be proposed. The lattice plays a distinguished role, and arises naturally in the 
ansatz for the discrete fermion field y/*^. We find: the discretisation satisfies the necessary algebraic 
relations, the reduction of the Dirac equation is earned out explicitly, and by choosing the ansatz 
i/A^, y/^, one finds the naive and staggered ys actions after application of a lattice projection 
operator 1^ to the "true action" then we discuss briefly the question of the chiral anomaly for 
the theory based on y. At the end, we discuss gauging and sketch the applicability to SUSY. 

1. Mathematical tools 

As announced, a central role is given to the Whitney map W . It associates a differential 
form to a chain, so that a vertex is mapped to a zero-form, an edge to a one form and a square 
to a two-form. In order to express the resulting forms, we need coordinates, these are given by 
the hypercubic generalization of the standard simplex. In each 2D elementary cell a, (2) of the 
complex K, we introduce local coordinates {x,y), in the interval [0, 1] x [0, 1]. For the present 
purpose, the main point about the Whitney forms is that given a face CJy,(p) of a given square CJ,- (2)> 
the Whitney form W(ay (2)) is defined over the entire coordinate region associated to <7,-,{2)- For 

' The details will be given in a forthcoming article 
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example^ VK([01]) = (1 —y)dx is defined over the entire interval [0,1] x [0,1]. This simple fact 
is conceptually important, the cells of highest dimensionality serve as local open sets in which 
operations can be defined. The discrete version of the exterior derivative is such that 

dW = Wd^, (1.1) 

which means that if the space of Whitney forms is considered, the discrete operation is exact, at 
every point in the interval [0, 1] x [0, 1] corresponding to the square CJ,- (2)- 

For the discrete Hodge star, the commonly used discrete operator is only approximate. A dual 
complex L is introduced, defined by joining the barycenters of edges and squares of the original 
complex K (image of K under -k^), of course, the resulting cells a^^^j overlap with four different 
squares cj^p)' there is no way the operation can be exact. Nevertheless, by using a dual lat- 
tice, the co-homology is correctly captured and the Poincare duaUty H^p\M) = H^^^p\M) indeed 
holds. A discrete wedge is introduced (N{p,q) is a normahzation) 

<^i,ip)^'^^J-il)=^(P^l)ll^k,{P+l) / W{0^p))AW{0j,(^^)). (1.2) 

k •'^k.[p+q) 

The discrete wedge is non-associative, for the simplicial version it is due to the presence of nu- 
merical factors, for the hyper-cubic complex, it is dependent on the chains having overlapping 
support, since if n / then CJj (^j a'^ ^i,(q) = in contrast Whitney forms always have 
overlapping support in the interval [0, 1] x [0, 1] corresponding to CJ,- (2)- 

We now supply this scheme with a contraction operation |^. It is defined using a Cartesian 
product on chains with a particular interpretation, it is 

j>a(,) =v^ X T]„/r(a(,)), (1.3) 

where both the vector and contracted chain T],,A:(a(,.)) are mapped to Whitney forms but the chain 
on the left of x is interpreted as a function coefficient. The Whitney map is then modified to 

W(ax/3) = 9(°)(lV(a))AW(j3). (1.4) 

The map deletes the form part of W(a). This product x is compatible with the other opera- 
tions. A key property is that the contraction is exact and the product is local in the sense that both 
W{o) and VK(j3) should have overlapping support to give a non-zero result. 

To start describing what is new in this talk, we propose to replace the discrete Hodge star 
which maps one complex to its dual by one which uses the product space x instead. Let 

*^a,;(p) = a,-(p)xa^('>), (1.5) 

where a'^^'^P^ maps under W to one of the globally constant forms^ {1, dx^ , dx^, dx^ Adx^}. The 
new discrete * is exact, and does not require a double complex. Previously, the presence of the 

^For the square [0123], [0] has coordinates (0,0), [1] is (1,0), [2] is (1, 1) and [3] is (0, 1). 
■'For example W{g^) = dx^ . 
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dual lattice L doubled the degrees of freedom, a form of species doubling^. Finally, the contraction 
operation may be used to define a discretized Clifford product (/ and g are chain degree dependent) 

^a.) vay,(,) =I(-l/((-l)^/|_ .../f^^a,(,)) A^/f^^ a,,^^^ (1.6) 
p 

The discrete wedge is approximate, but it is exact when one of the two chains maps to a constant 
form, this is the case with the chiral projection operators P^^\ 

2. Derivation of the lattice actions 

Let us start by writing down the natural action, 

^ = [ W{^'^)A*W{{d-5)(p^), (2.1) 
Jm 

having then specified the inner product, what remains is to find how to parametrize 0^. In the 
continuum theory, the field is an in-homogeneous form i.e 

0^ = + (j)idx^ + (l)2dx^ + (pudx^ A dx^. (2.2) 

By application of the chiral projection operator P'^^, the equation is reduced to the Dirac equation, 
the two-component spinor )//^(') = i/aVP'^Ms then 

1^(1) = (^>i (1 - idx^ A dx^) + ^2idx^ - idx^)- (2.3) 

The projection can be carried out in the present model as we will see, and the discrete DK operator 
is now exact i.e 

W{d'^ -5^) = {d-5)W. (2.4) 

The ansatz for the discrete fermion field is the key to deriving the lattice actions. We note that 
with the new discrete Hodge star, it is mandatory to use the product x for the basis of fields, since 
otherwise it can be shown that 5<7 = for any chain a. The ansatz then is the following discrete 
version of Eq. where the coefficients A(/ y) are placed on the lattice (labelled by a pair of 
integers (ij)), and 

V^"" = I H j) [('■ ' ^•)] - ( I [(', j)] ) X + Aj. y) [(/, j)] X (a • - ia^) , (2.5) 

we then verify it is invariant under P^'^ and substitute this ansatz into the DK equation 

(J^-5^)i/A^ = 0. (2.6) 

This leads to four equations, of the form i^k^k<yk) x a''. The chains on the left of x form two 
pairs of identical equations. The two independent equations are 

I(^(UlJ)-4j))[O^i)('>l'i)]-''(4;+l)-470)[(''-/')(''^ =0 (2.8) 

('■,;•) 

^also the new ★ playing the role of 'f projects correctly to the lattice under 3^, below. 
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These equations, after application of W , appear to be approximate in that they are functions of 
\,x^ and have no non-trivial solutions in the coefficients Nevertheless we can now connect 
to the lattice formulation, by noting that if we evaluate W{d^ — 5'^)y^ at a given vertex [(/, j)], we 
find precisely the naive 2D Dirac equation on the lattice 

a^5^Vfj)=0, (2.9) 

where is the symmetric difference operator. The lattice action can then be derived by inserting 
the projection operator £P = L(i,;) 8{^~X(ij)) in the action (x^^j j'j are the coordinates of i.e 

S^F= [ ^W{\if^)Ai.W{{d^-5^)Y^) (2.10) 
Jm 



where is the ansatz Eq.2^. For the staggered fermion, we modify the ansatz to Ys below. Of 
course one needs to introduce a double spacing for spinors, this may be done by introducing a dual 
lattice (same as L but with the new Hodge star) with its own collection of operations (A, J,*, /y)^, 
defined as we did in K. The ansatz is then (the Kronecker delta identify even and odd sites): 

V^l = I Sii2)flSj^2).o^j)[{iJ)]-i{I, 5,(2),i5,-(2).iA;^.)[(/,7)]) x a'^ (2.11) 

(IJ) {iJ) 
+ X {Si(2),lSj(2)flO^ - '5,-(2),o5y(2),ia^), (2.12) 

after subdividing the cells into the collection of finer ones B = KHL, one can write the action 
Ys) and after application of the lattice projection operator ^ one obtains the staggered 
fermion action. 

Let us now take a different viewpoint, we do not insert the projection operator ^ in the action 
=5^, and we evaluate ^(y/^, y^) to get a matrix M such that 

^(V/^, y"") = £ Y{ij)MY{k,i), (2.13) 
(U) 

it contains more couplings than the projected action namely diagonally opposite vertices in 

a square are coupled. This still respects the locality, and is chiral. On can add a Laplacian term 

Y/\*id5 + dd)Y (2.14) 

which breaks chirality and then project onto the lattice to a Wilson term. It is important to note 
that the lattice theory defined with the operator M does not have the same symmetries as the action 
^ which has continuous symmetries. This brings the possibility of using the action with the 
lattice parametrization of the field space y'^- In this case one can write explicitly the chiral current 
in differential geometric language. Since then chiral symmetry is not broken at the classical level, 
there is a strong suggestion that the anomaly may be derived between lattice and continuum in the 
framework presented here (there is no momentum cut-off). 

We now discuss gauging. The gauge field is a one-form, so it is associated to edges as 

= + l.i)] +Aii,MiJ+i)[i'J){iJ+y)], (2.15) 
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this makes the parallel with the link variable formulation easy to obtain. We consider the non- 
compact gauging 

^G= I W(v/^) A★W((^/*^-5^+A^V^)^/^^), (2.16) 

JM 

after inserting the projection operator, and identifying the leading terms in the expansion of the link 
variable U = e'^, we find (showing only ;c-axis link variables, k aie lattice unit vectors) 

(ij) (iJ) k 

the second term being not gauge invariant, we must impose a path ordering to prevent such term. 
It is then clear that the slightly non-local nature of the wedge product (within a cell) enables one 
to extract the link variable. Note that the latter cannot be used as a Wilson type term since it is not 
gauge invariant. An alternative non-compact treatment is to do the minimal coupling without using 
the wedge product but rather using the product x . Such approach also enables one to consider pure 
gauge theory and one finds that = F. 

We conclude with our prescription for lattice SUSY. A link between the DK fermion and ex- 
tended twisted SUSY (which is a topological theory) has been made [^. The coiTcsponding lattice 
theory is link variable based, and at an algebraic level is described by non-commutative geometry 
[0]. Our proposal also works at the algebraic level: the super-parameters are not discretized (e.g 
Wess-Zumino action), the partial derivatives are replaced with d^^ = (d^a) x a*', where d'^ is the 
jU-axis component of d^. The DK part of the twisted action is as defined above, including both 
chiral projections. This con^esponds to D = N = 2. For example. 

One immediately verifies that the algebra 

{Q, Q^} = id^, {<2, Q^} = -ie^vd\ (2.19) 
is realized on the field i//^ = Ys + where the flavor is identified as a twisted SUSY index. 
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